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ABSTRACT 

Context. B[e] supergiants are luminous, massive post-main sequence stars exhibiting non-spherical winds, forbidden lines, and hot 
dust in a disc-like structure. The physical properties of their rich and complex circumstellar environment (CSE) are not well under- 
stood, partly because these CSE cannot be easily resolved at the large distances found for B[e] supergiants (typically i 1 kpc). 
Aims. From mid-IR spectro-interferometric observations obtained with VLTI/MIDI we seek to resolve and study the CSE of the 
Galactic B[e] supergiant CPD-57° 2874. 

Methods. For a physical interpretation of the observables (visibilities and spectrum) we use our ray-tracing radiative transfer code 
(FRACS), which is optimised for thermal spectro-interferometric observations. 

Results. Thanks to the short computing time required by FRACS (< 10 s per monochromatic model), best-fit parameters and un- 
certainties for several physical quantities of CPD-57° 2874 were obtained, such as inner dust radius, relative flux contribution of the 
central source and of the dusty CSE, dust temperature profile, and disc inclination. 

Conclusions. The analysis of VLTI/MIDI data with FRACS allowed one of the first direct determinations of physical parameters of 
the dusty CSE of a B[e] supergiant based on interferometric data and using a full model-fitting approach. In a larger context, the study 
of B[e] supergiants is important for a deeper understanding of the complex structure and evolution of hot, massive stars. 



Key words. Stars: individual: CPD-57" 2874; Methods: observational, numerical; Techniques: high angular resolution, interferomet- 
ric; Stars: emission-line, Be 
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_^J. | 1. Introduction phase of ongoing or recently ended rapid mass transfer and dust 

; . formation. 
Q\ The B[e] stars form a heterogeneous group composed of ob- 

O jects at different evolutionary stages, but with many similar ob- One of the B[e] sub-classes is composed of luminous 

O servational signatures characterising the so-called "B[e] phe- (log(L*/L©) £ 4) post-main sequence objects: the B[e] super- 

^ ■ nomenon": simultaneous presence of hydrogen emission lines, giant stars (hereafter sgB[e]). Previous spectroscopic and po- 

j> \ low-excitation forbidden and permitted metallic lines in emis- larimetric observations of sgB[e] (e.g. Zickgraf et al. 1985; 

• *-j . sion, and a significant infrared (IR) excess mainly caused by hot Magalhaes 1992) show that the wind of these massive and lumi- 

/\ ■ circumstellar dust. Another common property of B[e] stars is the nous stars is composed of two distinct components: (1) a wind 

H ! presence of a non-spherical circumstellar environment (hereafter of low density and high velocity and (2) a wind of high den- 

■ - - ■ CSE; e.g. Zickgraf 2003). sity and low velocity. Zickgraf et al. (1985) proposed a picture 

Lamers et al. (1998) defined five B[e] sub-classes, one of where the S S B W winds consist of a hot and fast radiation-driven 

which contains unclassified stars. Miroshnichenko (2007) pro- w / nd and a s1 ™' mu J ch r i ool ^ r an 1 d 1 denser a fact ° r of 

posed an additional B[e] sub-class (the FS CMa stars) to explain 10 or 10 ) equatorial wind. This disc-like structure provides a 

at least part of the unclassified B[e]-type stars as binaries at a natural explanation for the existence of dust around those ob- 

jects, because the presence of dust requires regions of suffi- 
ciently high density and low kinetic temperatures. One possible 
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these discs (e.g. Meynet & Maeder 2006). Owing to their physi- 
cal characteristics (fast rotation, disc-like CSE, high luminosity, 
evolutionary status), it has also been suggested (Vink et al. 2009) 
that sgB [e] might share evolutionary links with rapidly rotating 
O-stars and long-duration gamma-ray bursts (GRBs). 

Because of the large distances of sgB[e] (<: 1 kpc for the 
closest ones), the geometry and physical structure (e.g. density 
and temperature distribution) of their CSE could be only quite 
recently directly probed, thanks to modern high angular resolu- 
tion (HAR) techniques. For example, Domiciano de Souza et al. 
(2008) used ESO's VLT/VISIR instrument to directly measure 
the typical size of the dusty CSE of the sgB[e] MWC 300 from 
diffraction-limited mid-IR images. Among HAR techniques, op- 
tical/IR long baseline interferometry (OLBI) provides the high- 
est resolving power, allowing the study of sgB[e] CSEs at angu- 
lar resolutions ~ 1 - 10 milliarcseconds (mas). For example, 
Millour et al. (2009) combined adaptive optics (VLT/NACO) 
and OLBI (VLTI/ AMBER and VLTI/MIDI) to detect a compan- 
ion around the sgB[e] HD 87643 and to perform an extensive 
study of the dusty CSE of the binary system. 

In the examples above as well as in most works based on 
HAR data of sgB[e], two different strategies are commonly 
adopted to interpret the observations: (1) geometrical analyti- 
cal modelling and (2) radiative transfer modelling (e.g. using the 
Monte Carlo approach). The geometrical analytical models have 
the advantage to be very fast to calculate, allowing a full model- 
fitting (for example a x 2 minimisation) and error estimate of the 
model parameters (mostly geometrical parameters of the CSE). 
However, these simple models do not give access to physical 
parameters of the targets such as temperature and density distri- 
butions, optical depths, etc. On the other hand, most radiative 
transfer models present a consistent description of the physi- 
cal conditions of the CSE. However, because these models are 
quite complex, they demand a lot of computing time, which pre- 
vents one from exploring a large domain of the parameter space 
and also from obtaining a good estimate of the uncertainties on 
the fitted parameters. In this work we adopt a third approach 
for the data modelling, which tries to keep the advantages of the 
other approaches, without the drawbacks. To this aim we use our 
fast ray-tracing algorithm for circumstellar structures (FRACS), 
which is based on a parametrised CSE combined to a simpli- 
fied radiative transfer (no scattering). A complete description of 
this algorithm is given in Niccolini, Bendjoya & Domiciano de 
Souza (2010; hereafter paper I). 

In the present paper we apply FRACS to study the CSE of the 
Galactic sgB[e] CPD-57° 2874 (also named Hen 3-394, WRAY 
15-535) based on mid-IR spectro-interferometric observations 
performed with ESO's VLTI/MIDI beam-combiner instrument. 
Previous near- and mid-IR interferometric observations of CPD- 
57° 2874 directly revealed an elongated CSE that is compatible 
with a disc-like structure formed by gas and dust (Domiciano 
de Souza et al. 2007; hereafter DS07). However, because only 
a limited number (four) of baselines was available and since 
the authors adopted simple analytical models, only geometri- 
cal parameters could be derived from this first analysis of CPD- 
57° 2874. As shown below, the use of FRACS allowed us to con- 
firm the previous results and, most importantly, to derive physi- 
cal parameters for this Galactic sgB[e]. 

In Sect. |2] we give the log of the VLTI/MIDI observations 
and describe the data reduction procedure. In Sect. [3] we pro- 
vide a new distance estimate of CPD-57° 2874 obtained from 
spectroscopic observations with FEROS. A short reminder of the 
ray-tracing code FRACS is presented in Sect. |H followed by the 
results obtained from a model-fitting analysis of the VLTI/MIDI 



observations (Sect.|5]l. A discussion of the results and the conclu- 
sions of this work are presented in Sects. |6]and|2l respectively. 

2. VLTI/MIDI observations 

The interferometric observations of CPD-57° 2874 were per- 
formed with MIDI, the mid-infrared 2-telescope beam-combiner 
instrument of ESO's VLTI (Leinert et al. 2004). All four 8.2 m 
unit telescopes (UTs) were used. The N-band spectrum as well 
as spectrally dispersed fringes have been recorded between A 
7.5 /im and A 13.5 fim with a spectral resolution of R ^ 30 us- 
ing a prism. In total, «b = 10 data sets have been obtained with 
projected baselines (B pro j) ranging from =* 40 m to =* 130 m, and 
baseline position angles (PA) between 8° and 105° (from 
North to East). A summary of the VLTI/MIDI observations of 
CPD-57° 2874 is given in Table [1] and the corresponding uv- 
plane coverage is shown in Fig. Q] 

Table 1. Summary of VLTI/MIDI observations of CPD- 
57° 2874: data set index, date, Coordinated Universal Time 
(UTC) of observation, baseline configuration, projected baseline 
length and position angle. 



# 


date 


^obs 


UT config. 




PA 






(UTC) 




(m) 


CO 


1 


2004-11-01 


08:51:55 


UT2-UT4 


85.1 


37.8 


2" 


2004-12-29 


05:52:12 


UT2-UT3 


45.2 


18.6 


3~ 


2004-12-29 


07:26:06 


UT2-UT3 


43.9 


35.1 


4° 


2004-12-31 


06:04:03 


UT3-UT4 


54.8 


79.6 


5" 


2004-12-31 


08:02:48 


UT3-UT4 


60.9 


104.8 


6 


2006-11-09 


07:15:36 


UT1-UT4 


129.7 


8.6 


7 


2006-12-13 


08:35:55 


UT1-UT3 


93.8 


28.6 


8 


2006-12-31 


08:13:05 


UT1-UT4 


125.8 


63.5 


9 


2006-12-31 


08:59:06 


UT1-UT4 


121.8 


72.5 


10 


2007-01-05 


08:35:09 


UT1-UT3 


88.1 


42.8 



" Data previously used by DS07. 




Fig. 1. uv-plane coverage: 
projected baselines (length 
and position angle) for the 
VLTI/MIDI observations of 
CPD-57°2874 (further de- 



180° 165 tails are given in TableQ}. 

The MIDI data were reduced with the MIA+EWS data re- 
duction package, which includes two different sub-packages: 
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Table 2. Observation log and angular diameters of calibrators 
(values from DS07) used to derive the calibrated N-band visibil- 
ities ofCPD-57°2874. 



date 


fobs 


UT coring. 


Calibrator 


©UD 




(UTC) 






(mas) 


2004-11-01 


09:28:16 


UT2-UT4 


HD94510 


2.16 ± 


0.11 


2004-12-29 


04:12:26 


UT2-UT3 


HD37160 


2.08 ± 0.20 


2004-12-29 


05:29:32 


UT2-UT3 


HD37160 


2.08 ± 


0.20 


2004-12-29 


06:13:08 


UT2-UT3 


HD50778 


3.95 ± 


0.22 


2004-12-29 


07:47:21 


UT2-UT3 


HD94510 


2.16 + 


0.11 


2004-12-31 


02:15:59 


UT3-UT4 


HD50778 


3.95 ± 


0.22 


2004-12-31 


03:04:33 


UT3-UT4 


HD50778 


3.95 ± 


0.22 


2004-12-31 


06:31:19 


UT3-UT4 


HD94510 


2.16 ± 


0.11 


2004-12-31 


07:19:17 


UT3-UT4 


HD 107446 


4.54 ± 


0.23 


2004-12-31 


07:41:22 


UT3-UT4 


HD94510 


2.16 ± 


0.11 


2006-11-09 


07:45:21 


UT1-UT4 


HD94510 


2.16 ± 


0.11 


2006-12-13 


02:04:34 


UT1-UT3 


HD23249 


2.33 ± 


0.01 


2006-12-13 


08:10:27 


UT1-UT3 


HD94510 


2.16 ± 


0.11 


2006-12-13 


08:54:41 


UT1-UT3 


HD94510 


2.16 ± 


0.11 


2006-12-31 


06:58:10 


UT1-UT4 


HD94510 


2.16 ± 


0.11 


2006-12-31 


07:44:55 


UT1-UT4 


HD94510 


2.16 ± 


0.11 


2006-12-31 


08:37:22 


UT1-UT4 


HD94510 


2.16 ± 


0.11 


2007-01-05 


08:08:04 


UT1-UT3 


HD94510 


2.16 ± 


0.11 



the MIA package developed at the Max-Planck-Institut fiir 
Astronomie, and the EWS package developed at the Leiden 
Observatory^]. 

While MIA is based on the power spectrum analysis, which 
measures the total power of observed fringes (Leinert et al. 
2004), EWS coherently adds the fringes after correction for opti- 
cal path differences (instrumental as well as atmospheric delays) 
in each scan (Jaffe et al. 2004). The data reduction results ob- 
tained with the MIA and EWS packages agree well within de 
uncertainties. 

The instrumental transfer function at each spectral chan- 
nel was obtained from the observations of calibrator stars with 
known uniform-disc diameters (®ud)- The calibrators used in 
the data reduction and the adopted angular diameters and uncer- 
tainties are listed in Table|2] The calibrated visibilities were cal- 
culated from the ratio of the targets' raw visibilities and the aver- 
age transfer function derived from the calibrator measurements 
of the corresponding night (Fig. 0. The error of the calibrated 
MIDI visibilities is of the order of 5% - 10% and includes the 
raw visibility error as well as the error of the transfer function. 
The uncertainties on the calibrator angular diameter are negligi- 
ble compared to the standard deviation of the transfer function. 
Usually, ^3-5 calibrator measurements per night were avail- 
able (Table 0. In the few cases where only one suitable cali- 
brator observation was available we assumed a typical transfer 
function error of 5% to estimate the errors on the final calibrated 
visibilities. 

VLTI/MIDI also provides spectral fluxes of CPD-57 2874 in 
the N-band. On average, all fluxes are compatible within 10% 
with those previously presented by DS07 (MIDI and ISO-SWS 
spectra). We here considered the VLTI/MIDI spectrum used by 
DS07 with an uncertainty of 20% (Fig. 0). This larger uncer- 
tainty ensures a complete agreement between the observed MIDI 
and ISO fluxes at all wavelengths. We note that the mid-IR spec- 
trum of CPD-57 2874 does not show any clear evidence of an 
important silicate feature around 10 fim. 



1 The MIA+EWS software package is avail- 
able at [hrtp://wwwmpia-hd.mpg.de/MIDISOFT/| and 
http://www.strw.leidenuniv.nl/" nevec/MIDI/index.html. 
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Table 3. Equivalent widths (EW) of the Na I absorption lines (AA 
5890 A, 5896 A), obtained from our FEROS data. The relative 
uncertainty of these measurements is about 20%. 



Line 5890A 5896A 

abs. comp. 1 abs. comp. 2 abs. comp. 1 abs. comp. 2 
W(A) OH O80 O05 O70 



3. FEROS observations and distance estimate 

In addition to our mid-IR interferometric observations, we 
obtained high-resolution optical spectra of CPD-57° 2874. 
The spectra were recorded with the high-resolution Fiber-fed 
Extended Range Optical Spectrograph (FEROS), attached to the 
2.2-m telescope at ESO in La Silla (Chile). FEROS is a bench- 
mounted Echelle spectrograph with fibers, covering a sky area 
of =* 2" x 2" and a wavelength range from 3600 A to 9200 A. 
Its spectral resolution is R 55 000 (around 6000 A). We have 
adopted its complete automatic online reduction, which includes 
the heliocentric correction. The FEROS spectra were obtained 
on 2008 December 21. We recorded two exposures of 1000 sec- 
onds with S/N of ~ 60 in the 5500 A region. 

These observations were used to estimate the distance of 
CPD-57° 2874. A previous distance estimate of d = 2.5 kpc has 
been proposed by McGregor et al. (1988), assuming that this star 
belongs to the Carina OB association. 

From our FEROS high-resolution spectra it is possible to es- 
timate the distance of CPD-57° 2874, based on the statistical re- 
lation cited by Allen (1973). This relation uses the equivalent 
widths of interstellar Nal lines. In our data, each Nal line is 
composed of two absorption components. However, owing to 
the lack of data from different epochs, it is impossible to see 
any temporal changes, which would allow us to derive a possi- 
ble circumstellar contamination. We have therefore assumed that 
both Nal components are of interstellar origin. The measured 
equivalent widths are given in Table [3] Our estimated distance 
for CPD-57° 2874 is d = 1.7 kpc with an uncertainty of 0.7 kpc. 
This large error is firstly due to a possible contamination from 
the circumstellar emission component and the saturation of the 
absorption one, and secondly to a systematic error caused by the 
statistical relation used. Within the error bars, our distance es- 
timate is roughly compatible with the result of McGregor et al. 
(1988). We have considered both distances in our analysis: 1.7 
and 2.5 kpc. 



4. Description of the model FRACS 

Here we present a short description of our numerical model 
FRACS and the parametrisation adopted to describe a sgB[e]. 
A full description of FRACS is given in paper I. 

FRACS is based on the ray-tracing technique using quadtree 
meshes for the grid and the full symmetries of the problem to 
hand to significantly decrease the computing time necessary to 
obtain monochromatic images (calculated with 300 x 300 pixels) 
within seconds. Complex visibilities and fluxes can be directly 
derived from these monochromatic images. FRACS neglects 
scattering, an approximation well suited to interpret spectro- 
interferometric observations in the thermal IR. Indeed, com- 
pared to absorption of light from dust, scattering can be ne- 
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glected in the IR and beyond (paper I; Lamers & Cassinelli 
1999). 

To analyse the VLTI/MIDI data of CPD-57°2874 we 
adopted the same parametrised description of a sgB[e] (central 
star and dusty CSE) as given in paper I. Below we summarise the 
main equations of this description in axis-symmetric spherical 
coordinates, and define the free parameters used in the model- 
fitting. 

We assume the specific intensity from the central regions of 
the star to be a power-law with spectral index a and level at 
a fiducial wavelength Aa = lOjum: 



(i) 



This emission from the central region includes a contribution 
from the stellar photosphere and from the continuum radiation 
(free-free and free-bound) of the circumstellar ionised gas. The 
spectral index a is sensitive to the nature of the central source. 
Panagia & Felli (1975) and Felli & Panagia (1981) give theoret- 
ical values for the spectral index for spherical envelopes: a =* 4 
for a blackbody and a 2.6 for a fully ionised gas (free-free 
emission) with an electron density proportional to r~ 2 . Their es- 
timates are valid within the Rayleigh-Jeans domain of the spec- 
trum, which fits to our case when we consider the hot central 
parts of a sgB[e] in the mid-IR. 

A radius R s = 54-Rq was adopted for the central region. This 
value is used simply as a scaling factor and to convert / s lo to the 
observed 10/im flux from the central region: 



(2) 



We note that at distances of a few kpc the central regions of a 
sgB[e] are not resolved by VLTI/MIDI and can thus be consid- 
ered as point sources. 

For the gas number density we adopt the bi-modal distribu- 
tion used by Carciofi, Miroshnichenko, & Bjorkman (2010) to 
study another B[e] star. Similar density distribution descriptions 
were adopted for Be stars (e.g. Stee et al. 1995). The adopted 
distribution is motivated by the two-wind scenario proposed by 
Zickgraf et al. (1985) and assumes a fast polar wind, a slow equa- 
torial outflow, and a latitude-dependent mass loss rate. The num- 
ber density of dust grains is therefore given by 



1 + A 2 \ 1+Ai (sine)" 



AW 1 + A 2 (sin0)" 



(3) 



where (r, 0) are the radial coordinate and co-latitude, and n m is 
the dust grain number density at 6 — 90° and at r = R m , which 
is the inner dust radius where dust starts to survive. Dust is con- 
fined between R m and the outer dust radius R out . The value of 
R out cannot be determined from the VLTI/MIDI data and has 
been fixed to a high value: 750 AU (the exact value does not 
affect our results). 

The parameter A i controls the ratio between the equatorial 
and polar mass loss rates. More precisely, (1 + A\) specifies the 
ratio between the mass loss rate per unit solid angle in the equa- 
tor and the pole. In the bi-modal scenario, the poles are assumed 
to be much less dense than the equator, therefore A\ must have 
a large value (k. 10). Our models indicate that for a wide range 
of values this parameter does not have a strong influence on the 
results because it can always be compensated by n m (see also dis- 
cussion in paper I). Therefore, we have arbitrarily fixed (1 + A\) 
to 50. 



The parameter A 2 indicates how much faster the polar wind 
is compared to the slow equatorial wind. (1 +A 2 ) is the equatorial 
to polar terminal velocity, i.e., v oo (90 o )/v oo (0°). This parameter is 
also quite uncertain and in principle can assume values ranging 
from ~ 1 to ~ 100. We kept A 2 as a free parameter, although it is 
not well constrained from the observations as shown in the next 
sections. 

Finally, parameter m controls how fast the mass loss (and 
consequently the density) drops from the equator to the pole. 
Defining the disc opening angle A#d as the latitudinal range 
within which the mass loss rate is higher then half its equato- 
rial value, we have 



A#d = 2 arccos 



2A, 



I 1 

2 arccos | — 



(4) 



High m values correspond to thinner regions of high den- 
sity around the equator. These dense, slowly flowing disc-like 
regions around the equatorial plane of sgB[e] stars provide 
favourable conditions for dust to form and survive. Different ap- 
proaches exist to define regions of dust formation (e.g. Carciofi, 
Miroshnichenko, & Bjorkman 2010). Here we adopt the rela- 
tively simple assumption where dust is allowed to exist only 
within the disc opening angle, i.e., at co-latitudes between 90° - 
O.5A0 d and 90° + 0.5A6» d . 

The dust grain opacity was calculated in the Mie theory (Mie 
1908) for silicate dust and for a dust size distribution follow- 
ing the commonly adopted MRN (Mathis, Rumpl & Nordsieck 
1977) power-law oc a -3 - 5 , where a is the dust grain radius. The 
Mie absorption cross sections are computed from the optical in- 
dices of astronomical silicate (Draine & Lee 1984; see also pa- 
per I). One possibility to reproduce the absence of a silicate fea- 
ture in the N-band spectrum of CPD-57 2874 is to have rela- 
tively large grain sizes. We thus used grain radii ranging from 
a = 0.5 to 50 //m, which can relatively well reproduce the ob- 
served spectrum (Fig. @}. We checked that ignoring scattering 
remains a valid assumption for this dust distribution with large 
grains. By neglecting the dust albedo the visibilities and fluxes 
are affected by only a few percent (;S 3.5%) within the N-band 
(further details in paper I). 

The temperature structure of the dusty CSE is given by 



T(r) = T m 



(5) 



where T m is the dust temperature at the disc inner radius R m , 
i.e., the dust sublimation temperature. To be consistent with our 
choice of dust composition we require that T m < 1500 K. The 
coefficient y is expected to assume values % 1 . 

Finally, because OLBI is sensitive to the projection of the 
object's intensity distribution onto the sky, there are two angles 
related to this projection: 

- the inclination of the disc plane towards the observer i (0° 
for pole-on view and 90° for equator-on view). 

- the position angle (from North to East) of the maximum 
elongation of the sky-projected disc PA d . This angle is de- 
fined for i + 0°. 

Thus, the 10 free parameters (nf ree ) of the model are: P., a, 
Tin, J, Rm, i, PAd, A 2 , n m , and m. 



5. Model-fitting with FRACS 

Here we use FRACS with the parametrised sgB[e] description 
defined in the last section in order to interpret the VLTI/MIDI 
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observations of CPD-57° 2874 through a model-fitting proce- 
dure. 

To ensure spectrally independent observations for the model- 
fitting we decided to consider one data point every 0.5 /im 
between 8 yum and 13 yum. This step approximately corresponds 
to twice the used spectral resolution width (AA = A/R). 
Additionally, to avoid poor visibility calibration owing to the 
Earth signature of ozone around 9.6 fim we have not included 
observations in this spectral region in our analysis. Finally, 
the same spectral sampling (n^ = 10 wavelengths points) 
was adopted for the VLTI/MIDI visibilities and spectrum. This 
choice also provides faster calculations because it is not neces- 
sary to compute model images at too many wavelengths. 

We have performed ax 2 minimisation simultaneously on the 
VLTI/MIDI visibilities and fluxes using a Levenberg-Marquardt 
(LM) algorithm (Markwardt 2008). In order to treat the visibili- 
ties and fluxes on the same level (similar weights) we have min- 
imised a x 1 like quantity defined as (see further details in paper 
I): 

where V° bs and Vk are the observed and modelled visibility mod- 
ulus for baseline index k, F° bs and F ; are the observed and mod- 

j 1 

elled mid-IR fluxes for wavelength index j. <Tv,k and cr F j are the 
estimated errors on the visibilities and fluxes. 

The starting parameter values for the fit were determined 
from physical considerations of the CSE and from the previous 
results from DS07. Below we consider the reduced x 2 defined 
by x 2 = X 2 /(2 n B n A - "free), where «f ree = 10. The LM algorithm 
stops when the relative decrease in^ is less then 10~ 3 . For the 
CPD-57 2874 data, the LM algorithm reaches the^ minimum 
(xtn-n r ) m a f ew hours (=; 2 - 3 h) on a single CPU. 

Figure [2] shows the intensity map of the model correspond- 
ing to x\( m r (best-fit model) for our distance estimate of 1 .7 kpc, 
which also corresponds to the lowest x 2 ■ The visibilities and 

1 A min,r 

fluxes for the best-fit model are shown, together with the obser- 
vations, in Figs. [3] and |4] These plots show that the model well 
reproduces most observations within their uncertainties for both 
adopted distances (1.7 and 2.5 kpc). In particular the slightly 
curved shape of the visibilities is well reproduced by FRACS. 
The models indicate that this curved shape is probably caused 
by the combined fact (1) that the intensity maps have different 
relative contributions from the central source and from the dusty 
CSE at different wavelengths, (2) that the optical properties of 
the adopted dust grains are wavelength-dependent even if there 
is no strong silicate feature seen in the spectrum, and (3) that 
the angular resolution significantly changes along the observed 
wavelengths. 

The model parameters at x\av, r an d their uncertainties are 
listed in Table |4] The derived parameters are almost indepen- 
dent of the adopted distance, except of course for those scaling 
with the distance. The uncertainties of the parameters have been 
estimated from^ maps calculated with 21x21 points for each 
pair of free parameters (45 pairs) and centred on the x 2 nin r posi- 
tion. All 45 Xt ma P s ar e shown in Figs. [5]to[7]for d = 1.7kpc. 
These maps show that the^ space presents a well defined^ in , 
without showing several local minima in the parameter domain 
explored. Additionally, they provide visual and direct informa- 
tion on the behaviour of the model parameters in the vicinity 




Fig. 2. Intensity map of CPD-57° 2874 at 10/um for the best-fit 
FRACS model obtained for a distance d — 1.7 kpc (see Table|4]i. 
The image scale is in log of the specific intensity R. 

of x 2 ■ > revealing, for instance, potential correlations between 
certain parameters. 

We have estimated the parameter uncertainties in a conser- 
vative way by searching for the maximum parameter extension 
in all^j 2 maps corresponding to;^^. + Ax 2 , where Ax 2 = 1 (see 

contours in Figs. [5] to [7j- This choice of Ax 2 sets a lower limit 
confidence region of =s 60% to the parameter uncertainties. This 
limit results from two extreme assumptions about the data: 

- Data points per baseline are completely dependent (corre- 
lated): because the same set of stars is used to calibrate all 
visibilities of a given baseline, we can consider a limiting 
case where all these visibilities are correlated. This assump- 
tion implies that only 10 independent visibility observations 
are available (this corresponds to the number of baselines). 
The flux at each spectral channel can still be considered to be 
independent. This pessimistic assumption leads to the lower 
limit of =* 60% to the formal confidence level for Ax 2 = 1, 
corresponding to only 20 independent observations (10 base- 
lines and 10 fluxes). 

- All data points are completely independent (uncorrelated): 
an upper limit of 100% of formal confidence level is ob- 
tained if we assume that all data points are independent. 
Then the uncertainties derived from Ax 2 = 1 are very con- 
servative (overestimated). 

Hence the parameter uncertainties given in Table [4] corre- 
spond to a confidence level of at least 60%, but most probably 
they are somewhat overestimated. 

In the next section we present a physically motivated discus- 
sion of the derived model parameters of CPD-57 2874. 

6. Discussion 

6.1. Geometrical parameters (PA d , R in , and i) 

Let us first compare the derived geometrical parameters (PAd, 
Ri n , and i) with those previously obtained by DS07 from ellipti- 
cal Gaussian models fitted on a sub-set of the VLTI/MIDI data 
used here. 

The geometrical parameter PAj can be directly compared 
with the major-axis position angle of the ellipse previously deter- 
mined by DS07 (^ 143° - 145°). As expected, the two estimates 
of PAd are identical within their error bars. 
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Fig. 3. VLTI/MIDI visibilities of CPD-57°2874 (circles) and the best-fit FRACS visibilities obtained from a x 1 minimisation for 
a distance d = 1.7 kpc (solid curve) and d = 2.5 kpc (dashed curve). The model-fitting was performed simultaneously on the 
visibilities and spectral flux. The visibilities effectively used for the fit are shown as filled circles together with the corresponding 
visibility error bars. The good quality of this fit is reflected by a reduced x\m °f ~ 0-55 for both distances (see details in Tabled. 



Because the bulk of the thermal IR emission comes from 
the internal regions of the disc, one can expect the inner dust 
angular radius (R^/d) to be comparable to (or slightly smaller 
than) the major-axis half width at half maximum (HWHM) of an 
elliptical Gaussian. Indeed, the Ri n /d derived here agrees with 
the major-axis HWHM (= 0.5FWHM) given by DS07: 4.5 < 
HWHM (mas) < 8.0. 

Domiciano de Souza et al. (2007) estimated a CSE view- 
ing angle i ~ 30° - 60° from the minor- to major-axis ratio of 
the fitted elliptical Gaussian model. This estimate agrees fairly 
well with the more precise determination of this parameter given 
here. 

The comparison of these parameters shows on one hand that 
the geometrical parameters obtained with FRACS agree with 
those obtained from a simpler approach using analytical models. 
On the other hand, this comparison clearly shows that, at the cost 
of somewhat higher but similar computing times, FRACS gives 
us access to physical parameters of CPD-57°2874 that cannot 
be extracted from simple geometrical analytical models. 



6.2. Continuum emission from the central source (II and a) 

In contrast to what may be initially expected, we show below 
that the central source (star and continuum emission components 
such as free-free, free-bound) contributes to almost half of the 
total mid-IR radiation of CPD-57° 2874. The total lOyum flux of 
the best model for d = 1 .7 kpc is F tot = 7.9 x 10~ 13 W irT 2 yizrrT 1 
(see Fig. The contribution of the dust CSE alone to the 10 /mi 
flux computed with FRACS is F A = 4.4 x 10 _13 WrrT 2 /an -1 . It 
follows therefore that the 10pm flux from the CSE is =* 56%F tot 
and from the central source is 44%F (ot . Similar results are 
obtained for d = 2.5 kpc. 

We note that in the particular case of CPD-57 2874 where 
i =i 60° and where the dust is confined in a relatively narrow 
disc (opening angle ~ 7°), these relative flux contributions can 
be directly obtained from F\ (derived from I\; see Table |4j>. 
This is valid if there is no absorption of the central regions by 
the CSE. 
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Table 4. Best-fit model parameters and uncertainties derived for 
CPD-57° 2874 from a x 2 minimisation. The uncertainties were 
estimated from the^ 2 maps. 



Adopted distance 


d = 


1.7 kpc 


d = 2.5 kpc 


Reduced X 2 min 




= 0.54 


x 2 - = 

/x i"nin r 


0.56 


Model parameters 


value 


error 


value 


error 


I s , (10 5 Wirr 2 um-' str 1 ) 


2.2 


+0.7 
-0.7 


4.2 


+ 1.8 
-1.4 


a 


2.4 


+ 1.3 
— 1.2 


2.4 


+1.4 
-1.4 




1498 


+1042 


1500 


+ 1050 
—535 


y 


1.02 


+0.71 

—0 29 


0.86 


+0.43 
—0 24 


R m (AU) 


12.7 


+3.6 
-2.9 


14.4 


+5.1 

-4 


iC) 


61.3 


+ 10.8 
-18.2 


59.6 


+ 11.8 
-21 2 


PA d O 


140.3 


+ 12.3 
-14.0 


139.4 


+ 13.7 
-15.1 


Ail 


-0.98 




-0.98 




n m (m- 3 ) b 


0.30 




0.33 




m b 


332 




377 




Other derived parameters 


value 


error 


value 


error 


F s Aa (1(T 13 WnT 2 yum-')^ 


3.5 


+i.i 
-i.i 


6.8 


+2.8 
-2.3 


R,„/d (mas) 


7.5 


+2.1 
-1.7 


5.8 


+2.0 
-1.6 


A0 d (°) 


7.5 




7.0 





" To be physically consistent with the adopted dust composition the 
upper limit for 7" in is 1500 K, even though the^ 2 maps were allowed to 
explore higher temperature values. 

b Not well constrained. 

c Observed mid-IR flux from the central region at Aq = 10 fim (Eq.[2). 



14 - 




8 9 10 11 12 13 

wavelength Qurn) 

Fig. 4. VLTI/MIDI flux of CPD-57°2874 (thick solid grey 
curve) and the +20% adopted uncertainty (dots). The thin solid 
and dashed curves are the best-fit model fluxes for assumed dis- 
tances of d = 1.7 kpc and d = 2.5 kpc, respectively. The wave- 
lengths used for the fit are those from Fig. [3] 



Although the uncertainties on the spectral index a are rela- 
tively high (^ 50%), the derived value (= 2.4) suggests an im- 
portant contribution of free-free continuum radiation from the 
central regions (Felli & Panagia 1981). 

6.3. Temperature structure of the dusty CSE (T m and y) 

The dust temperature at the inner radius T m is found to be 
=* 1500 K (the imposed upper limit for the fit), which is consis- 



tent with the definition of T m itself and with the chosen silicate 
dust composition. The large upper limit uncertainties in T m and 
the fact that the best-fit Ti n is 1500 K for d = 2.5 kpc indicate that 
a Tin slightly higher than 1500 K could still be compatible with 
the observations. A higher T m value is consistent with, for in- 
stance, different dust compositions. However, we prefer to keep 
our choice of dust composition and to have a relatively large up- 
per limit uncertainty in T^, since the present observations do not 
provide strong constraints on the exact dust composition. 

The derived value for the coefficient of the temperature pro- 
file y indicates an almost linear decrease of the disc temperature 
as a function of r. The steepness of this temperature profile lies 
between those expected for a non irradiated, adiabatically cool- 
ing disc (T °c r~ 4 ^ 3 ) and a reprocessing disc (T oc r~ 3 ^ 4 ) (c.f. 
Porter 2003). This implies that a non negligible part of the repro- 
cessed radiation from the inner parts of the disc escapes without 
being re-absorbed, so that the disc cools down faster than highly 
optically thick discs, without of course reaching the limit of a 
purely adiabatic cooling. 

6.4. Parameters related to the density law (n in , A 2 , and m) 

Our results suggest that the observed mid-IR visibilities and 
fluxes cannot strongly constrain each individual parameter re- 
lated to the density law: n m , A%, and m. The x 2 maps show a 
significant correlation between these three parameters, indicat- 
ing that they are degenerated for the available VLTI/MIDI ob- 
servations. Even if their uncertainties are significant, the values 
of these parameters at / ^ 2 nin suggest 

- a low inner density n m corresponding to relatively low CSE 
optical depth in the mid-IR (^ 0.2 along the line of sight 
around lO/zm). 

- a polar-to-equatorial terminal velocity ratio 
Voo(0°)/v co (90°) = 1/(1 + A 2 ) - 50, compatible with 
the values found in the literature (e.g. Zickgraf 2003). 

- a high value for m, translating into a quite narrow opening 
angle (< 10°) for the dust disc. 

Because n; n , A 2 , and m are not well constrained, we fitted the 
observations by fixing these parameters to their values in Table|4] 
in order to investigate their influence on the remaining parame- 
ters. We have also fixed the J 1 ;,, to 1500 K. The fit was performed 
for d - 1.1 kpc, starting from slightly different values from those 
in Table [4] . The x 1 an d values obtained for the free parameters 
(I\ , a, y, R m , i, PAd) are essentially the same as in Table|4](dif- 
ferences are only a small fraction of the parameter uncertainties). 
We have also checked that the uncertainties on the other param- 
eters are not affected by the fact that n; n , A 2 , and m are not well 
constrained. 

6.5. Data analysis from a model with fewer free parameters 

Thanks to the data analysis performed here we found that some 
parameters of the sgB[e] model adopted for CPD-57° 2874 can- 
not be well constrained from the available VLTI/MIDI data. Of 
course this is not necessarily a general conclusion because it 
depends on the nature of the studied target and on the spectro- 
interferometric data available. 

The results from our analysis of CPD-57° 2874 indicate that 
it is justified to consider a simplified version of the model de- 
scribed in Sect [4] As shown in the previous section, the param- 
eters related to the density law are the less constrained by the 
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Table 5. Best-fit model parameters derived for CPD-57° 2874 
from a x 2 minimisation using a model with fewer free parame- 
ters then the initial one (c.f. Sect. 16.51 ). Uncertainties were esti- 
mated from the Levenberg-Marquardt (LM) algorithm and can 
be considered as lower limits to the errors on the derived param- 
eters (see discussion in the text). 



Adopted distance 
Reduced^ 


d = 1.7 kpc 
y 2 = 54 


t/ = 2.5 kpc 
xLns = 0-56 


Model pcircimeters 


value error 
LM 


Vcilue error 
LM 


I\ (10 5 Wm- 2 yimr 1 str 1 ) 
a 
7 

R m (AU) 
iC) 

PA d O 
n m (itT 3 ) 

A0 d O 


2 1 +al 

-0.1 

7 4 +0.2 
92 +0 07 

-0.07 

1 1 +2 

ll.U _ 20 

60.5 +}j 
139.8 +|;° 
09 +0 06 

-0.06 
7 c +4.4 
'■ J -4.4 


4 +o.i 

^ -0.1 

7 4 +0.2 

0.85 tS 
13.9 ! 2 * 
59.3 +{j 
139.3 

1 1 +0 07 

1 1 -0.07 
t q +4.4 
-4.4 



data. Let us thus consider an alternative density law where the 
number density of dust grains is given by 

n(r g)= i n m (^) 2 ; 90° - O.5A0 d < 9 < 90° + O.5A0 d (J 
\0 ; 9 < 90° - 0.5 A9 d and 9 > 90° + 0.5 A6» d 

The parameters A i , A2, and m are not present in this simpler 
density prescription. Only n m and A#d are necessary to define 
the density structure. Based on our previous results we have also 
fixed T m to 1500 K. The number of free parameters is thus re- 
duced from 10 to 8, namely, I\ , a, y, R\ n , i, PAd, n m , and A#d. 

Following the previous procedure we have performed a x 2 
minimisation on the visibilities and fluxes using a LM algorithm. 
The results are shown in Table [5] Note that the best-fit param- 
eters completely agree (well within the uncertainties) with the 
previous values obtained with the initial model (Table|4]l. 

Only n m shows an important difference compared to the pre- 
viously tested model. Contrarily to the previous density law 
(Eq. 0, the new one (Eq. [7} assumes that the density is con- 
stant along a given r inside the dust disc. In order to obtain sim- 
ilar fluxes and optical depths as before, n m has to be somewhat 
smaller than the previous value. 

Table [5] also provides the parameter uncertainties estimated 
with the LM algorithm, which are smaller than those estimated 
from the^ maps (given in Table|4]i. These smaller errors appear 
because all data points are assumed to be independent and no 
covariance matrix is used in the error estimations. The estimated 
errors from LM in Table[5]can then be considered as lower limits, 
while the errors from the x 2 maps can be considered as upper 
limits to the parameter uncertainties (see discussion in Sect. [5). 

7. Conclusions 

The dusty CSE of the Galactic sgB[e] CPD-57° 2874 was spa- 
tially resolved thanks to mid-IR spectro-interferometric observa- 
tions performed with the VLTI/MIDI instrument. Several physi- 
cal parameters and corresponding uncertainties of this star were 
derived from a x 2 minimisation and from the analysis of x 2 
maps. The physical quantities derived include the inner dust ra- 
dius, relative flux contribution of the central source and of the 



dusty CSE, dust temperature profile, and disc inclination (refer 
to Tableland Sect.0for details). 

To our knowledge, this is the first direct determination of 
physical parameters of the dusty CSE of a B[e] supergiant based 
on interferometric data and using a model-fitting approach from 
ax 2 minimisation. This was possible thanks to FRACS, which 
adopts a parametrised description of the central object and of 
the dusty CSE combined to a simplified radiative transfer (no 
scattering). Its main advantage is computation speed (< 10 s per 
monochromatic image with 300 x 300 pixels). Because it is fast, 
FRACS allows us (1) to explore a large parameter space domain 
for a true global x 2 minimisation and (2) to more realistically 
estimate the uncertainties on the best-fit parameters. We would 
like to recall that contrarily to a model such as FRACS, simple 
geometrical models do not allow a simple and direct access to 
physical parameters and uncertainties of the dusty CSE. 

Future complementary observations could be included to 
measure new CSE parameters and/or to reduce the uncertainties 
and the correlations on some parameters that were not strongly 
constrained by the VLTI/MIDI observations alone. Consistently 
with the domain of validity of FRACS, these complementary 
data should be obtained at wavelengths above the near-IR, where 
the dust starts to contribute to a significant amount of the stellar 
flux (UV to visible observations cannot be consistently modelled 
by FRACS). For example, the closure phase information from 
the MATISSE beam combiner (second generation instrument on 
VLTI; Lopez et al. 2009, Wolf et al. 2009) could inform on the 
disc inclination and opening angle, which are parameters that 
can influence the symmetry of the mid-IR intensity distribution. 
Also, near-future interferometric observations in the millimetre 
with ALMA will be more sensitive than MIDI and MATISSE to 
the thermal emission from the colder regions of the dusty CSE. 
Thus, these observations will probably better constrain the tem- 
perature profile of the dust and also provide direct information 
on the structure and actual size of the CSE, allowing also for 
mass estimates. Moreover, this information on the colder dust is 
very important for the study of the evolutionary history of sgB[e] 
stars (e.g. mass and angular momentum losses, chemical enrich- 
ment and interaction with the close interstellar medium). 
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Fig. 5. Reduced^ 2 maps from VLTI/MIDI observations of CPD-57° 2874 for all 45 combinations of free parameters of our model 
(FRACS). The maps are centered on the^ iinr (= 0.54) position and correspond to the distance of 1.7 kpc. Contours are drawn for 

X 1 + A^ 2 , with Ax 2 = 0.3 ,1,3. The parameters are ordered as in Table|4] The map scale is given in logarithmic units. Details on 
the model-fitting procedure are given in Sect. [5] 
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Fig. 6. Continuation of Fig. [5] 




Fig. 7. Continuation of Fig. [6] 



